Control of Surface Wettability via Strain Engineering 
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Abstract 

(— I . Reversible control of surface wettability has wide applications in lab-on-chip systems, tunable 

^ I optical lenses, and microfluidic tools. Using a graphene sheet as a sample material and molecular 

dynamic (MD) simulations, we demonstrate that strain engineering can serve as an effective way 
to control the surface wettability. The contact angles 9 of water droplets on a graphene vary from 
72.5° to 106° under biaxial strains ranging from —10% to 10% that are applied on the graphene 
layer. For an intrinsic hydrophilic surface (at zero strain), the variation of 9 upon the applied 
strains is more sensitive, i.e., from 0° to 74.8°. Overall the cosines of the contact angles exhibit 
a linear relation with respect to the strains. In light of the inherent dependence of the contact 
angle on liquid-solid interfacial energy, we develop an analytic model to show the cos as a linear 
function of the adsorption energy Eads of a single water molecule over the substrate surface. This 
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'nT . model agrees with our MD results very well. Together with the linear dependence of Eads on biaxial 

^^ ' strains, we can thus understand the effect of strains on the surface wettability. Thanks to the ease of 

reversibly applying mechanical strains in micro/nano-electromechanical systems (MEMS/NEMS), 
• 1-^ 

rS ' we believe that strain engineering can be a promising means to achieve the reversibly control of 

- - - surface wettability. 
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I. INTRODUCTION 



Understanding the surface characteristics and controlhng the wettabihty of sohd surfaces 
are fundamental topics in chemical physics and serves as a basis in many applications [l|. The 
Lotus-Leaf-like super-hydrophobic surface, for example, has attracted enormous attention 
in the past decades because of the superior water-repellent and self-cleaning properties |2|. 
On the other hand, super-hydrophilicity of surfaces has also seen itself wide applications, 
e.g., improving the filtration efficiency of polymer filter thin films |3| and the boiling heat 
transfer efficiency in heat pipes |j, |5|. Recently active and reversible control of the surface 
wettability is becoming a very attractive research topic and various applications are proposed 
and demonstrated. 

Engineering the solid surface structures is the most widely used method to control the 
wettability, e.g., via introduction of surface roughness/defects and changes of the chemical 
properties |2|, y, Ig, |7| . It is reported that the surface energy of graphene increased as more 
defects were induced, leading to a hydrophilic nature |8|. Rafiee et al. found that the 
contact angle of graphene film can be controlled when the graphene film was dispersed 



on a substrate by performing hig 
proportion of acetone and water 



i-power ultra-sonication in the solvent with a controlled 
9| . Despite the great success of this approach, there are 
several trade-offs. Inducing roughness, defects and chemical groups will affect the structural 
integrity and thus reduce the robustness of the performance in practice (e.g., under harsh 
environments) [10|, together with some other undesired consequences, e.g., enhanced slip 
friction force for the motion of fiuid droplets, and so on. Additionally, most of the methods 
to change the solid surface structure and wettability are usually permanent and lack of active 
controllability and reversibility, which are highly desired in novel micro/nano-fiuidic devices 
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Reversible controls of wettabilit y h ave been recently demonstrated by some elegant meth- 
ods, including the light irradiation 12|,|l3|, the electrochemical surface modifications |l4l.ll5|. 
applying an electric field 161, and the thermal treatment 171. But it still suffers some limita- 

nn 

tions, e.g., a small controllable range of contact angles (about 11° [3|, Il3[), a large hysteresis 
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■ Il6|. delayed dynamic motion of liquid droplets 12[, and a low number of life cycles 



13|. Most of these limitations can be attributed to the fact that the wettability changes are 



accomplished by the conformational transition of the surface molecular structures. 



In this paper, we will demonstrate that the strain engineering can serve as an effective 
way to reversely control the wettability of an atomically smooth surface without damaging 
the structure. Graphene is selected as the solid substrate, because of its superior mechanical, 
electronic, and bio-compatib ility properties, which render it an ideal material building-block 
in the nano-fluidic devices 18|. We studied the wetting behaviour of water droplets on 
graphene sheets with a biaxial strain e range from —10% to 10% using molecular dynamics 
(MD) simulations. The contact angles 9 of water droplets on the graphene surface vary from 
72.5° to 106° in this strain range. For the intrinsic hydrophilic surfaces (i.e., at zero strain), 
the variation of 9 upon the applied strains is even more sensitive than that of intrinsic 
hydrophobic surfaces (e.g., from 0° to 74.8°). Overall we find a linear relationship between 
the cos 9 and e. In the end, an analytical model will be presented to explain this observation. 

II. SIMULATIONS AND METHODOLOGY 

Our molecular system is illustrated in Fig. [U^a): a water droplet on a single graphene 
sheet. The graphene sheet is biaxially strained, ranging from —10% to 10%. Accordingly 
the carbon-carbon bond length ace changes from 0.9acCo to l.lacco! where acco denotes 
the carbon-carbon bond length of strain-free graphene (e.g., acco = 0.142 nm). We placed 
four different sizes of water droplets on the strained graphene sheets, i.e., 748, 1885, 3709, 
6600 water molecules, respectively. In our MD simulations, the positions of carbon atoms 
are fixed. We have tried flexible graphene substrates with applied 0% and 10% biaxial 
strain and found a less than 5° difference in the calculated contact angles from those of the 
corresponding rigid substrates. 

All MD simulations were performed with the LAMMPS code [19J. A time step of 1.0 
fs was used and the total simulations time was about a few nanoseconds. We used the 



CHARMM force field and the SPC/E model |20| for water with the SHAKE algorithm 
2l|. The van der Waals (vdw) interactions between water molecules and carbon atoms 
were described by a Lennard- Jones (LJ) potential between oxygen and carbon atoms, i.e. 
0(r) = 4e((cr/r)^^ — (cr/r)^). The vdw forces were truncated at 1.2 nm with long-range 
Columbic interactions computed using the particle-particle particle-mesh (PPPM) algorithm 
|22| . Water molecules were kept at a constant temperature of 300 K using the Nose-Hoover 
thermostat. 
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FIG. 1. (Color online) (a) A water droplet on strained graphene where the carbon atoms are fixed; 
(b) Density map of a water droplet with 6600 water molecules on graphene sheet at 0% strain, 
which has been averaged along the radial direction; (c) Determination of the contact angle 9 and 
the base curvature l/rs by fitting a circle to the free surface of the water droplet in (b). 



While maintain the same hexagonal lattice with graphene for the monolayer solid sub- 
strates, we select three sets of LJ parameter for the water-solid interactions, corresponding 
to macroscopic contact angles of 91.2°, 52.1°, and 133° on the zero strained solid surfaces. 
There are two motivations. First, we can study the strain engineering effects on the con- 
tact angles for substrates with different wettability. Second, the contact angle of graphite 
measured in experiments is scattered, varying from 0° to over 115°. It should be noted 
that value of about 90° is commonly accepted 23|, u3\- The parameters of a = 0.3190 nm 



and e = 4.063 meV were adopted to reproduce such an angle (231]. We kept the same a 



value but altered the parameter e as 5.848 meV or 1.949 meV to represent the hydrophilic 
or hydrophobic strain- free solid substrates, respectively 23|. We label these substrates as 
graphene, hydrophilic-surface and hydrophobic-surface, accordingly. 

It usually took a few hundred picoseconds to reach equilibrium and the simulations were 
then continued for two more nanoseconds to collect data. From the MD molecular trajecto- 



ries, water density maps were obtained by introducing a 3-diniensional grid with the size of 
each cell as 0.5A x 0.5A x 0.5A. Reducing the 3D mesh into a 2-dimensional density map by 
averaging along radial direction leads to Fig. dJ^b), which is the density map of the droplet 
with 6600 water molecules on a graphene sheet with 0% strain. We further averaged the 
2D density maps of 10000 frames of the trajectory in a total of 2 ns duration. To obtain 
the water contact angle from such a map, a two-step procedure was adopted following the 



reference 23||. First, the boundary of the droplet surface was determined within every single 
layer that was vertical to z direction by using 0.5 g/cm as a critical density 25|. Second, a 



circular best fit through these points was extrapolated to the solid surface where the contact 
angle Q was measured as shown in Fig. [U^c). Note that the points of the surface below a 
height of 8 A from the solid surface were not taken into account for the fit, to avoid the 
influence from density fluctuations at the liquid-solid interface. 

III. RESULTS AND DISCUSSIONS 

Fig. |2] summarized the contact angles of the water droplets with different sizes on 
the strained sheets of graphene (Fig. [2]J^a-b)), hydrophilic-surface (Fig. [2t^c-d)), and 
hydrophobic-surface (Fig. HJ^e-f)). The measured microscopic contact angles Q depend 



on the size of the droplets. According to the modified Young's equation 23|, l26|, the 



macroscopic contact angle ^oo can be related to the microscopic contact angles as: 

T 1 
COS^ = COS^oo (l) 

where r is the line tension, l/r^ is the fitted base curvature of the water droplets, and 
7^^ is the interfacial free energy per unit area for the liquid-vapor interface. In Fig. [21 the 
cosines of contact angles measured in our MD simulations follow the linear relation very well 
with respect to l/r^ (Eq.([l])). Extrapolating the linear relations in Fig. [2]|^a,c,e) yields the 
macroscopic contact angles Boo- 

Fig. El^b), (d), and (f) depict the contact angles d and 6*00 (hexagon symbols) as a 
function of the applied strains. The cosines of the (macroscopic) contact angles exhibit a 
linear relation with the mechanical strains. For the three types of surfaces, positive strains 
(stretching) always result in a decrease of cos Q and thus the increase of contact angle ^, and 
vice versa. The magnitude of the angle changes is, however, quite different for the three 
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FIG. 2. (color online) Cosine of the contact angle as a function of the base curvature l/rs of 
the droplets on (a) graphene, (c) hydrophilic-surface and (e) hydrophobic-surface; Cosine of the 
contact angle as a function of the biaxial strains applied on (b) graphene, (d) hydrophilic-surface, 
(f) hydrophobic-surface. 6 



types of solid surfaces. The 6'oo of graphene varies from 72.5° to 106° for the apphed strain 
from —10% to 10%. On the hydrophobic-surface, the 6*00 varies from 121° to 140°. For the 
hydrophihc-surface, the change is much more significant, from 0° to 74.8°. In comparison, 
other methods have a smaller controllable range of contact angles. For example, by using an 
electric field, the (receding) contact angles of water droplets on a low density monolayer can 
be reversibly controlled between 20° and 50° 16|. Subject to a light irradiation, the largest 



contact angle change on langmuir-b 



azobenzene polymer is about 11° 13 1. 



odgett films with photoresponsive fluorine- containing 
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FIG. 3. (color online) line tension fitted from Fig. [2] using Eq. ([T]) as a function of the applied 
strain on the substrate. 



Fitting the contact angle vs. droplet size to Eq. ([T]) yields the values of line tension r 
as well. The results are shown in Fig. |3l They are on the order of (1.0 — 3.0) x 10~^^J/m, 
depending on the wetting properties of the substrates and applied strains (Fig. [3]). The 
obtained magnitude of r is consistent with the theoretical predictions and some recent 
experiments 27|]. It is known that a substrate with a hig her surface atom density and better 
wettability usually has a larger line tension r 28|, |29|. Stretching the graphene and the 



hydrophobic-surface will reduce the surface atom density and the wettability, thus leading 



to a dropped line tension value as observed in Fig. [31 In contrast, Fig. |3] suggests that the 
hydrophilic-surface (a nearly constant line tension, 2.7 x lO^^^J/m) is almost independent 
on these two factors. It is worth noting that there is still a debate on the quantitative value 
of r and data reported often have differences in several orders of magnitude 25|, |26 |. 

In light of the inherent relation between the contact angle and the water-solid interfacial 
energy ''ysL-, in the foUowings, we will report a linear relation between the cos6'oo and the 
adsorption energy Eads of a single water molecule over a solid surface, which will help us 
understand the effects of strain on the cos^oo (Fig- E])- The Young-Dupre equation [30 1 
correlated the macroscopic contact angle and the work of adhesion: 



7iy(l + cos^oo) = VrMo 



(2) 



where W = {'Jlv + Isv — 1sl)Ao is the work of adhesion, Aq is the contact area, •ysv, Isl 
and 7iy are the interfacial free energy per unit area for the solid-vapor, solid-liquid, and 
liquid-vapor interfaces, respectively. 

In our molecular model, the water-solid interaction energy is zero when they are separated 
infinitely far away from each other. Thus we can calculate the work of adhesion W by 



summarizing the molecu^ 



surface atoms 



i8i,y, 



ar pair-wise interactions between the water molecule and the solid 



3l|, i.e., 



W = - 



E{h)p{h)A{h)dh 



(3) 



where E{h) is the interaction energy of a single water molecule with a distance h from the 
substrate, p{h) are number density of the water molecules per unit area, and A{h) is the cross 
section area of the water droplet. As usual, the summation of pair- wise vdw interactions 
between one water molecule and all carbon atoms on the monolayer surface can be replaced 
by an integration assuming the continuous distribution of carbon atoms 30| . 
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where Eads is the lowest point of the potential well E{h), which occurs at h = a: 

Sire , a ,o Svre , a ,r, Sne , a 



E, 



ads 



E{h = a) 



{-—y{l-2s) (5) 



5^3 flee bVS (l+£:)acco 5a/3 ctcco 

The value of Eads represents the adsorption energy of single water molecule over a solid 

surface, which is linearly correlated with the applied mechanical strain e on the solid surface 

in Eq. IBj. 
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Our previous MD study [32j showed that there is a depletion layer at the interface with a 

thickness of 0.2 nm, between 0.2 nm and 0.5 nm is the so called 'first water layer' with a peak 

density p about 2-3 times higher than the bulk value, between 0.5 nm and 1.0 nm is a second 

layer with density exhibit small oscillation around bulk water density, beyond 1.0 nm is the 

bulk water with a constant density po- This density profile p{h) has a weak dependence 

on the strains applied on the surface. A similar density profiles are also observed for the 

hydrophilic and hydrophobic surfaces in this study. Thus, we can perform the integration 

[Eq. ([3])] in two separate parts: 

1 /"'*" n /"°° 

l + cos0oo = / p{h)E{h)dh ^/ E{h)A{h)dh 

ILV J hi iLvAo Jhu 

p{h)E{h)dh - Atail 



iLV J hi 



iLVJh 



hu 



p{h) 



ha^ 2a 



10 



dh \ Bads - Atail (6) 



3/i4 3/iio_ 

where the hi and hu represent the lower bound of the first water layer and the upper bound 
of the second, respectively. 

In Eq. ([6]), the first water layer has a monolayer thickness so that its cross section 
area is approximately the droplet contact area Aq. The rapid decay of E{h) for h > a 
implies that the second term on the right side of Eq. ([6]) is much smaller than the first one. 
Indeed, our calculations showed that Atau is less than 5% of the first term. Because the a 
is a constant (LJ potential between water and carbon atoms), the liquid- vapor interfacial 
'Jlv is independent of the solid surfaces, and the change of water density p(/i) upon strain 
engineering has a small effect in the integral term, which we will discuss in detail later, we 
can conclude that cosine of the macroscopic contact angle 6*00 should exhibit an approximate 
linear relation with respect to the adsorption energy Eads- 

Fig. m summarizes the cos 6*00 directly obtained from our MD simulations (Fig. [2]) as a 
function of Eads for the differently strained graphene, hydrophilic-surface and hydrophobic- 



surface. Indeed, it is a linear function, as predicted in Eq 



The fitted value of the slope 



is — 22.0/eV. The contact angle data from Werder etc. 23|] are also included in Fig. [Hand 
the agreement is good. Since Eads is a linear function of the applied strain [Eq. ([5])], we can 
understand the linear relation between cos 6*00 and e depicted in Fig. O 

It is interesting to investigate the changes of water structures upon the strains of graphene 
layer and how such changes would determine the contact angle. In our previous publication 
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FIG. 4. (color online) Cosine of the macroscopic contact angle 9oo as a function of the adsorption 
energy of a single water molecule on top of the strained graphene, hydrophilic-surface, hydrophobic- 



surface. We also included the simulated microscopic contact angle results of Werder etc. 



M- 



32], we analyzed the water structure (density, radial density function and structure factor) 
at the interface of the strained 'real' graphene. The two-dimensional radial distribution 
and structure factors of the first water layers on the —10%, 0% and 10% strained graphene 
substrates are almost identical, while the water density profiles along the perpendicular 
direction of the substrates are similar except a ±15% change in the peak density of the first 
water layers. We have carried out analysis on MD simulations results on hydrophobic and 
hydrophilic surfaces. Similar conclusions are obtained. We found that the change of water 
density upon strain engineering has a small effect in the integral term in Eq. (Q. For 'real' 
graphene, the integral term varies from — 22.2/eV to — 20.7/eV at a strain from —10% to 
10%. For hydrophobic surface, it varies from — 18.9/eV to — 17.0/eV. And for hydrophilic 
surface, it is from — 22.5/eV to — 22.2/eV. All of them are close to the fitted slope — 22.0/eV 
in Fig. m This shows our analytical model is consistent quite well with our MD results. This 
also indicates the change of Eads upon strain engineering plays a dominant role on cosine of 
the contact angle. 



10 



It is worth noting that in Ref. 23|, Werder et.aL correlated the contact angle 9^ to 
the equilibrium adsorption energy, in which they concluded that only in a certain range, 
the contact angle follows a linear relationship with the Eads- We believe that Eq. ^ (Fig. 
S]) is a better model to understand the relation between the surface wettability and the 
chemical/physical interactions of liquid and solid surfaces. 

Since our model is derived without presumption on the atomic structures of the substrate, 
it can be easily applied to 3-dimensional substrates (e.g. graphite) by summing the interac- 
tion energy between water molecules and the different atomic layers of substrate, although 
the surface atomic density of a substrate and the averaged Eads might not necessary be 
linear functions of the applied strain/stress. It is worth noting that since the Lennard- Jones 
interaction decays rapidly over distances, the interfacial energy is mainly determined by 
the first atomic layer of substrates (9|, |33| . We believe our results can provide some valid 
indications for a broad class of three dimensional substrates. 



Strain up to 30% can be readily exerted on graphene in experiments 34l-l36| . For example, 
an epitaxial strain of ~ ±1% builds up in the graphene when it is grown on different 
I; uniaxial strain up to ~ 1.3% to a graphene monolayer can be applied by 



substrates 



36 



using two- and four-point bending setups 35|; and uniaxial strain ranging from to 30% 



can be achieved for large-scale graphene films transferred to a pre-stretched substrate 
There are also many ways to reversibly control the strains applied on substrate materials 
by mechanical stresses or electric fields. For example, reversible strain obtained in Fe-Pd 
single crystals by com pre ssive stress-induced martensite variant rearrangement is reported to 
exhibit as high as 5% 40| ; Zhang etc. obtained a large reversible electric-field- induced strain 
of over 5% in BiFe03 films 41|. We, therefore, believe strain engineering can be a promising 
way to reversely control the surface wettability in practice. In practice, compressive strain 
often leads to ripples or folding of graphene. In this manuscript, the study on the compressive 
strain serves a theoretical interest. We aim to obtain a more complete picture of the relation 
between contact angle and molecular adsorption energy and substrate strain. 

Strain engineering, in our previous study, has been predicted to significantly change 
the slip length of water over a graphene layer 32|. It is thus interesting to investigate 
the correlation between the wettability and the slip length. It is often believed that a 
hydrophobic surface often has a larger slip length than that of a hydrophilic surface 42 |. 



Strain engineering the wettability and the slip length turn to be an exceptional case. In 
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Fig. |2l with a strain from —10% to 10%, the graphene is becoming more hydrophobic 
(e.g., contact angle increasing from 72.5° to 106°). The shp length obtained in our previous 
MD simulations was reduced from 175 nm to 25 nm. Liquid water can slip on a hydrophilic 
surface, even with a larger slip length than that on a hydrophobic surface 43|. To understand 
such a counter-intuitive observation, we should note that Eq. (E]) implies the dependence of 
contact angle on Eadsi which is an average of the van der Waals potentials over the whole 
solid surface (at h = a), whereas it is the energy barrier (i.e., the corrugation of the vdw 
potential profile) experienced by water molecules over a solid surface who defines the slip 



length 32|. Stretching a graphene layer results in a smaller E^ds [Eq. (E])] and thus a higher 



contact angle 6oq, but it leads to an enhanced energy barrier and thus a reduced slip length 



32| . This insight may help understanding the observed controversial correlations between 



the contact angle and the slip length in experiments 42 



44|. 



IV. CONCLUSIONS 

To summarize, using graphene sheet as a sample material, our MD simulations show 
that the wettability of a solid surface can be controlled by mechanical strains. Overall, the 
cosines of the contact angles exhibit a linear relation with respect to applied strains. For a 
graphene surface, the contact angles can be tuned from 72.5° to 106° under biaxial strains 
ranging from —10% to 10%. For droplets on intrinsic hydrophilic surfaces (at zero strain), 
the variation of the contact angle is more sensitive than that of a hydrophobic surface. 
To understand the strain engineering effect, we developed an analytical model to reveal 
a linear relationship between the cos6'oo and the adsorption energy Fads of a single water 
molecule over the substrate surface. The applied mechanical strains change the Eads and 
consequently alter the contact angle. We believe that the linear relationship between the 
cosine of the contact angle and the adsorption energy Eads is a general model to describe 
the surface wettability and the chemical/physical interactions of liquids and solid surfaces. 
Thanks to the ease of reversibly applying mechanical strains in MEMS/NEMS, we propose 
that strain engineering can serve as an effective means to achieve the reversibly control of 
surface wettability. 
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